Metaferroelectrics: Artificial Ferroelectricity in Metamaterials 
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Metamaterials are artificial periodic structures which represent effective homogeneous medium for 
electromagnetic fields. Here I show that there exists an important class of such composite systems, 
metaferroelectrics. They are properly designed metamaterials with non-ferroelectric constituents 
which, nevertheless, exhibit ferroelectric response. Furthermore, for relatively small samples the 
effect is anisotropic and strongly depends on the sample shape. Metaferroelectrics will help to 
understand deeper the physics of ferroelectricity and show enormous potential for applications in 
the field of nano-optoelectronics. 

PACS numbers: 77.80.-e, 77.84.Lf, 77.90.+k 



Metamaterials represent artificial periodic structures 
usually made of metallic nanoparticles. They exhibit a 
number of striking properties which lead to such unique 
effects as perfect lensing, electromagnetic cloaking, etc. 
[]]]. The lattice constants of these systems are smaller 
than the wavelength of electromagnetic waves. There- 
fore, for the electromagnetic fields, metamaterials effec- 
tively represent a homogeneous medium. Perhaps this is 
the main reason why the role of periodicity was not fully 
appreciated for tailoring their electromagnetic properties. 
In reality, as is demonstrated below, the proper choice 
of lattice parameters is of paramount importance, and in 
certain cases, it determines qualitatively new phenomena 
such as artificial ferroelectricity. 

Ferroelectrics have been a subject of intensive research 
for many years 0, B 0]- Nevertheless, the underlying 
physics of these materials remain less understood com- 
pared to their magnetic counterparts, ferromagnets [Bj]. 
More than a half century ago Slater proposed that depo- 
larization fields, stemming from the dipole-dipole inter- 
action, may be responsible for appearance of spontaneous 
electric polarization Q. However, it occurs that for ex- 
isting ferroelectrics dipolar fields alone do not suffice 
Here, I show that in metamaterials with non-ferroelectric 
constituents the depolarization fields lead to the sponta- 
neous polarization, and therefore, make possible to real- 
ize artificial ferroelectricity. Such metaferroelectrics will 
advance our understanding of physics of ferroelectricity 
and show a huge potential for possible applications 

In particular, let us consider metamaterial which con- 
sists of cubic lattice of metallic nanospheres of radius r 
embedded in a dielectric matrix with dielectric constant 
£ m . Clearly, the lattice constant d must obey d > 2r. 
To be specific the metal is assumed to be silver. Below 
the plasma frequency the dispersive properties of silver 
is best represented by the Drude dielectric function Q: 



of free electrons n t , as follows: 
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where the plasma frequency lo p is related to the density 
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and S is the damping constant. Further, (— e) is the 
electron charge, m cit is the effective electron mass in the 
metal, and eo is the permittivity of free space. 

The induced electric dipole moment p of an isolated 
silver sphere in the external electric field E is 
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This result is obtained for the static uniform electric field. 
However, it is still valid in the quasistatic limit when 
electromagnetic waves have wavelength much bigger than 
the sphere radius Q . Indeed, the comparison with exact 
Mie theory shows that for 3nm < r < 25nm the qua- 
sistatic approximation gives very accurate results. For 
the spheres with bigger radius quadrupole and higher or- 
der multipoles become important (loj . 

The dipole resonance frequency w R , at which p would 
become infinite in the absence of losses, is Re{e Dr (o; R ) + 
2e m } = 0. This gives: 
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It should be noted that, for r > 3nm silver nanoparti- 
cles, the values of w R obtained on the bases of equation 
(|4]) agrees very well with experimental data as well as 
with quantum mechanical calculations 0]. In the case of 
smaller spheres quantum size effects come into play and 
equation fl| is no longer valid. 

The dynamical e quat ion which correctly reproduces 
expression ([3]) reads [III , [ijj : 
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with p = — gx q is the absolute value of the oscillator 
charge q — 7V c e, and iV is the total number of free elec- 
trons in the sphere N e ~ (Airr 3 /3)n c . Further, x is the 
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average electron displacement. The damping constant 7 
and the mass of the oscillator m are some functions of fre- 



quency However, for our purpose only the values of 
these quantities in the static limit are important. At zero 
frequency 7 is an irrelevant parameter, and m = gN a m a{{ 
with g = ^/(3e m a;|). 

Now, in an artificial lattice the local field acting on 
a sphere is E + (L/eq)P, where L is the depolarization 
factor [7] , and P is the polarization density. The equation 
of motion ([5|) then becomes: 
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Taking into account that the unit cell volume is d 3 , the 
polarization density can be written as P = (—qx)/d 3 . 
This allows to rewrite equation ([6]) as: 
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with the effective frequency w off given by 
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Assuming that E ~ exp (iu>t) equation ([7]) gives: 
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The factor in front of E diverges at zero frequency when 
uj 2 it = 0. This means that arbitrarily small electric fields 
can induce nonzero electric dipole moment and the sys- 
tem becomes ferroelectric The condition that the 
metamaterial is ferroelectric, u>^ e = 0, can be cast in a 
more convenient form: 
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It is rather interesting that the lattice constant d F docs 
not depend on the parameters of the Drude dispersion re- 
lation |T]). Let us assume that the background dielectric 
is silicon e m = 11.68 (l3j |. For the (infinite) cubic lattice 
L = 1/3 [7], and equation (JTUJ) gives d F = 3.66r. Here the 
following remark is in order. If the metallic nanospheres 
were point objects, and in the case of highly symmetric 
cubic lattice, the dipole approximation employed above 
would be correct for arbitrary d F . However, the real 
spheres have a finite radius and the lattice constant must 
be bigger than certain d cr . For d < d CI coupling between 
neighboring spheres excites quadrupole and other higher 
order modes and the dipole approximation becomes inac- 
curate. The detailed numerical studies (see e.g. |l4ll5l|) 
show that d cr = 3.2r. Since d F > d al for the exemplary 
system considered here, the approximations involved are 
self-consistent and accurate. 

Moreover, as was stressed above, the spheres behave 
as effective dipoles. The dipole fields decrease with the 



distance as ~ 1/R 3 . On the other hand, the number of 
dipoles located at the distance R from the given lattice 
site can be estimated as ~ R 2 . Therefore, if the sample 
size is ~ R s , the field stemming from the surfice dipoles 
is proportional to ~ 1/R S , which is a slowly decreasing 
function. This means that, for relatively small systems, 
the depolarization fields depend on the sample shape [f| . 
In general, they are inhomogeneous and anisotropic, and 
therefore, are represented as a space-dependent tensor. 

Nevertheless, for ellipsoids, the depolarization fields 
appear to be homogeneous. If the coordinate axes co- 
incide with the principle axes of the ellipsoid, all non- 
diagonal elements of the depolarization tensor are iden- 
tically zero. Furthermore, the following normalization 
relation L x + L y + L z = 1 holds, where L x , L y , and L z 
are the diagonal elements of the depolarization tensor 0] . 

In the case of prolate ellipsoid, 
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with the major and minor axes a and b respectively (a > 
b), we have: 
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L y — L z — -(1 - Lx), 
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where the eccentricity e PE is given by e PE = y/l — (b/a) 2 . 
For a oblate ellipsoid, 



x 2 + y 2 z 2 



with the major and minor axes a and c respectively (a > 
c), the depolarization factors are: 



Lx — L y — -(1 - L z ), 
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L z = o- 22 (e OE - arctane OE ) , 
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and the eccentricity e OE is e OE = \J (a/c) 2 — 1. 
Then, the equation of motion |Q takes the form: 
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where the effective frequency tensor cD off reads: 
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Note that x and E in equation (|13|) are represented by 
the column vectors. 

The system exhibits the ferroelectric response in a 
given direction if the corresponding diagonal element in 
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equation (fl4|) vanishes. Below some important particular 
cases are considered in more details: 

(i) If a sample has the spherical shape L x = L y = L z = 
1/3. Therefore, in this case the depolarization fields are 
the same as for the infinite lattice and g? f p = 3.66r. The 
system response is isotropic and polarization density P 
is a three-dimensional vector P — (P x , P yi P z ). 

(ii) If a sample has the cylindrical or needle geometry 
which is oriented parallel to the x-axis L x = and L y — 
L z = 1/2. In this limiting case equations (fT0|) and (fl4|) 
give dp 1 = 4.19r. The system exhibits the easy-plane 
anisotropy and the polarization density is a planar vector 
P = (0,P y ,P z ). The ferroelectric response is isotropic in 
the plane normal to the x-axis. 

(iii) Finally, for a thin film or disk with the plane nor- 
mal to the z-axis L x = L y = and L z = 1. Thus, 
equations (10]) and HU) suggest that d® = 5.27r. The 
system shows easy-axis anisotropy, and therefore, the 
polarization density becomes a one-dimensional vector 
P = (0, 0, P z ) with only one nonzero component. 

In summary, here I have demonstrated that arti- 
ficial ferroelectricity can be achieved in properly de- 
signed metamaterials. Most interestingly, the con- 
stituents of the composite system are assumed to be non- 
ferroelectric. In particular, the suggested structure con- 
sists of a cubic lattice of silver nanospheres embedded 
in a silicon matrix. Moreover, it should be stressed that 
for relatively small systems, due to the long-range na- 
ture of dipole-dipole interaction, the corresponding lat- 
tice constant d F strongly depends on the sample shape 
and the effect becomes anisotropic. Metaferroelectrics 
show enormous potential for applications in the field of 
nano-optoelectronics . 
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